On the quantum mechanical nature in liquid NMR quantum computing 
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The quantum nature of bulk ensemble NMR quantum computing — the center of recent heated 
debate, is addressed. Concepts of the mixed state and entanglement are examined, and the data in 
a 2 qubit liquid NMR quantum computation are analyzed. It is pointed out that the key problem 
in the current debate is the understanding of entanglement in a mixed state system. The following 
points are concluded in this Letter: l)Density matrix describes the "state" of an average particle 
in an ensemble. It can not describe the state of an individual particle in an ensemble in detail; 2) 
Entanglement is a property of the wave function of a quantum particle(such as an molecule in a liquid 
NMR sample). Separability of the density matrix can not be used to measure the entanglement of 
mixed ensemble; 3)The evolution of states in bulk-ensemble NMR quantum computation is quantum 
mechanical; 4) The coefficient before the effective pure state density matrix, e, is an measure of the 
simultaneity of the molecules in an ensemble. It reflects the intensity of the NMR signal and has no 
significance in quantifying the entanglement in the bulk ensemble NMR system. We conclude that 
the liquid NMR quantum computation is genuine, not just classical simulations. 



Computers based on quantum mechanical principle of superposition can do computation much faster than classical 
computers jl]-^| ■ The proposals for quantum computing using liquid-state NMR have accelerated the experimental 
studies of quantum computing in demonstrating quantum algorithms and quantum error correction codes [^| p!s[ . 
However this approach has been in the center of a heated debate |l^ . The center of the debate is whether liquid 
NMR quantum computations carried out so far are genuine quantum computation. 

To implement a quantum computation, one has to use a quantum mechanical system to represent the wave function 
and to enbody the quantum computation operation. In a liquid NMR system at room temperature, the system is in 
a thermal equilibrium and the molecules arc not in the same quantum state. One has to construct an effective pure 
state by certain method, exhaustive averaging J25|] for example. After this, the density matrix for the system can be 
written as 

p=(l-e)I/d + e Pl , (1) 

where / is the identity matrix in d = 2 N dimension, and N is the number of qubit. p\ is the density matrix for a pure 
state. It is found that |]l9|-p2|| when e in equation (|l|) is smaller than a critical value which is true for all the present 
NMR experiments carried so far, the density matrix can be decomposed into a linear combination of products of the 
states of the individual qubits. Separable density matrix represents separable state which can be written as product of 
the composite particle wave functions and they are unentangled. It was argued |l9|-|2^] that since the density matrix is 
separable for small e, there was no entanglement. The bulk ensemble NMR computation experiments carried so far are 
not genuine quantum computations and can only be considered as classical simulations of the quantum computation, 



a point strongly disagreed by Lafiamme |23|. In this Letter we will address this issue by examining the concepts of 
mixed state and the concept of entanglement. We also analyzed the evolution of the density matrices in a 2 qubit 
NMR quantum computation and compare them with prediction of quantum mechanics. 

Microscopic particle, such as a molecule, obeys the laws of quantum mechanics. Its state is fully described by a 
wave function whose time evolution is the Schodinger equation. We do not talk about the state of a particle pure or 
mixed. Because it is always described by a wave function in quantum mechanics |24|| . Pure or mixed state refers to 
an ensemble. By pure state, we mean that all the molecules in an ensemble have the same quantum wave function. 
In other words, every particle in the ensemble can be described by the same wave function. By mixed state, we mean 
that the molecules in an ensemble are not in the same quantum state, but rather has a classical distribution of the 
molecules in different quantum states: n\ particles in state ri2 particiles in state \ip2) ■■■■ Each of the particle 
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in the ensemble is in a quantum state, described by a wave function. However, we do not have the detailed 
information for each particle and what we know is just the property of an average particle in the ensemble. To describe 
the property of an average particle in an ensemble, we use the density matrix which was invented by Von Neunmann 
p6| . In the 2 qubit system such as phosphoric acid for example, there are 4 quantum states, 1 00) , |01),|10) and |11), 
where the first bit represent the nuclear spin of the 31 P and the second bit represents the nuclear spin of hydrogen 
atom. If all the phosphoric acid molecules are in the same state, say, 1 00) , then we have a pure state. For liquid NMR 
sample at room temperature, the molecules are in the thermal equilibrium and the ensemble is described by a mixed 
state. 



The idea of entanglement was first formulated by Einstein, Podolsky and Rosen |27|. When we say the wave function 
of a quantum system is entangled we mean that the wave function of the whole system can not be factorized into a 
product of the wave functions of the constituent particles. It should be stressed that here the constituent particles are 
all quantum mechanically related one another, such as the nuclear spins of the hydrogen atom and 31 P in a phosphoric 
acid molecule, not the molecules of an ensemble. 

Now we study the concept of entanglement in a mixed state. First it is easy to see that the same density matrix can 
be prepared in numerous ways. For example, in a two qubit system, a density matrix which is 1/4 of a unit matrix 
can be prepared by putting into the ensemble with a quarter of particles in each of the 4 states with spin up and 
down along the z-axis: 

p = ||00)<00| + ~|01)(01| + j|10)(10| + ~|11)(11|. (2) 

The same density matrix can be obtained by putting a quarter number of particles into the ensemble in each of the 
4 Bell states. 

p = i(|00) + |ll))((00| + (ll|) + i(|00)-|ll))((00|-(U|) 

+5d 01 > + |10))((01| + (10|) + |(]01) - |10))((01| - (10|). (3) 

Though the density matrices are the same, the physical systems are completely different. There are infinite number 
of ways to prepare a system with a given density matrix. This is not surprising, since density matrix describes 
only the average property of a system. In this example, both ensembles can be viewed as an ensemble consisting 
of 4 sub-ensembles with each sub-ensemble in a pure quantum state. The molecules in these ensembles have quite 
different entanglement properties. In (^|), each individual molecule wave function is factorisable, whereas in (^), each 
individual molecule is non-separable and is enatngled. Since density matrix for a mixed ensemble describes only the 
average property and different physical systems can have the same density matrix, one can not tell whether a mixed 
state ensemble is entangled or not. Even if a density matrix is separable, it does not mean neccesarily the 
molecules in the ensemble is entangled unless the system is prepared in the way as suggested by the 
decomposition. Two ensembles such as (||) and (^) may have identical density matrices, but the physical systems are 
completely different, though we do not have any measurement to distinguish them. From the mathematical properties 
of the density matrix alone, one can not draw conclusions about the detailed properties of the individual particles 
of the system. The physical properties of the ensemble depend on the history or the method of its preparation p8| . 
The pure state is an exception, because all the molecules are in the same quantum state, the average property of the 
system is identical to the property of an individual particle |^9| . 

Secondly, even we know the history of the ensemble and locate the physical system, the density matrix of an 
ensemble reflect only the properties of an average particle in an ensemble and does not tell the information of an 
individual particle in the ensemble. From the density matrix, we can not say more about a particle than its average 
property. Taking for example the systems represented by (|^) and (||)(Here we use the expansion in the equation to 
specify its preparation procedure.) The difference in the properties of the invidual particles in the two ensembles can 
not be reflected in the density matrix. 

Thirdly, entanglement is a concept for a quantum system only, a single molecule for instance. We can talk about 
the entanglement of wave functions of nuclear spins of the hydrogen atom and the 31 P atom of a single phosphoric 
molecule. But we can not talk about the entanglement properties between two molecules in an ensemble. This is 
true for a pure ensemble. What is meant by entanglement is the entanglement within each individual molecule, not 
the entanglement between different molecules in the ensemble, and the molecules in the ensemble are 
independent of one another. The purity of the state is a measure of purity of the ensemble, or the synchronization 
or simultaneity of the wave functions of the molecules in the ensemble. When quantum computation is performed 
in a pure cnscmblc(a pure state), each molecule in the ensemble undergoes the same evolution simultaneously. For a 
mixed ensemble, the quantum computation operation is performed on molecules in different quantum states, including 
the state we arc interested in. Effective pure state preparation procedure is only a way to extract the required signal. 
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In a simple example, let's look at a one qubit system 
r ' () 

In a extremely simplified illustration, we can imagine that there are 



P = 







are in state 



Suppose the density matrix of the system can be written as 

2-level particle. Five of them 



0) and three particles are in state |1). When we measure the system which is a radio frequency pulse in 
NMR, both |0) — * 1 1) (upward) and |1) — > |0) (downward) transitions are induced. The transition signals meet in the 
coil, 3 upward transitions cancel with 3 downward transitions. The net signal in the coil is 2 upward transitions. This 
is equivalent to a pure system with two particles in state |0). When the number of qubits becomes large, effective pure 
state procedures are used to make unwanted signals cancel and the required signals collected. From this simple picture 
we see e, the coefficient before the effective pure state density in ([!]) is just a measure of the number of molecules in 
the required state whose signal can be read out. It is closely related to the intensity of the readout signal. There is 
no essential difference between an NMR experiment with e = 1 and one with a smaller e value. 

In a liquid NMR experiment, the evolution of the molecules in the ensemble are quantum mechanical. We have 
analyzed the experimental data between two steps in a phase matching study of the generalized quantum searching 
algorithm ||3(J. The details of the quantum algorithm is not important here. We give here the transformation matrix 
for a step in the quantum computation is 
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What we are interested in is the evolution of the density matrix, 
state tomography technique |3lf| directly from experiment, 



(4) 



The density matrix is obtained by the quantum 



P(l) = 



0.1794 0.1591 + 0.0208/ 0.0601 - 0.0001/ 

0.1591 - 0.0208/ 0.2453 0.1247- 0.0281/ 

0.0601 + 0.0001/ 0.1247 + 0.0281/ 0.3616 

-0.0483 + 0.0549/ -0.0514+ 0.1534/ 0.0099 - 0.0682/ 



-0.0483 - 0.0549/ 
-0.0514- 0.1534/ 
0.0099 + 0.0682/ 
0.2137 



(5) 



After one step of quantum computation manipulation, the density matrix is constructed by quantum state tomography 
again, 



0.2278 0.0858 + 0.0186/ 

0.0858- 0.0186/ 0.1006 

0.0640-0.0387/ 0.1019 + 0.0062/ 

0.0691 + 0.0372/ 0.1650 + 0.0893/ 



0.0640 + 0.0387/ 
0.1019- 0.0062/ 
0.3921 
0.0454 + 0.0111/ 



0.0691 - 0.0372/ 
0.1650- 0.0893/ 
0.0454-0.0111/ 
0.2794 



(6) 



The theoretical prediction of the density matrix using quantum mechanics through p' th = c ■ p ■ 

0.0891 



is 



Pth 



0.0891 
0.0758 - 
0.1146- 



0.1849 
0.0599/ 
0.0225/ 
0.0439/ 



0.0650- 
0.1377- 



0.0599/ 
0.0999 
0.0446/ 
0.0861/ 



0.0758 
0.0650 

0.0018 



0.0225/ 
0.0446/ 
0.3876 
0.0083/ 



0.1146- 0.0439/ 
0.1377- 0.0861/ 
0.0018 - 0.0083/ 
0.3277 



(7) 



The theoretical prediction of the evolution of the density matrix is in good agreement with the experimental mea- 
surement. 

We summarize the main points of this Letter: 1) The density matrix is not unique in locating a physical system. 
To uniquely specify a physical ensemble, one has to include the history of the ensemble. Separability of the density 
matrix when e is small alone does not mean the molecules in the ensemble are untangled. The molecules are entangled 
only when the ensemble is prepared in a way suggested by the decomposition. 2) Entanglement is a property of an 
individual particle in an ensemble, not a property of the ensemble as a whole. 3) The evolution of the density matrix 
in a liquid NMR system is quantum mechanical. 4) The procedures for preparing an effective pure state is only a 
way to extract the proper signals from the ensemble. The coefficient e of the effective pure state does not determine 
the entanglement property of an ensemble. It is only a measure of the intensity of the signals of the required state. 
This conclusion has an important practical implication for NMR quantum computation, a small e does not change the 
quantum nature of an NMR quantum computation. One can repeat the quantum computation for sufficient number 
of times to effectively enhance e. This will help in scale an NMR quantum computer. One can increase the number 
of repetition to compensate for the decrease of signals for large qubit numbers. 5) Finally we'd like to stress that 
power of the quantum computation lie both in the coherence of states (superposed states) and the ability to perform 
quantum mechanical unitary operation. Entangled state is a special case of the superposed state. During a quantum 
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computation, the quantum register is always in a superposed state, and sometimes it is entangled and sometimes it 
is unent angled. 

To conclude, the present NMR quantum computation is genuinly quantum mechanical, much more than a simula- 
tions of a quantum computation. Liquid NMR technique is still a very useful tool in experimental studies of quantum 
computation. 
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